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Abstract —Silicon Photomultipliers (SiPMs) are attractive can¬ 
didates for light detectors for next generation liquid xenon 
double-beta decay experiments, like nEXO (next Enriched Xenon 
Observatory). In this paper we discuss the requirements that 
the SiPMs must satisfy in order to be suitable for nEXO and 
similar experiments, describe the two test setups operated by the 
nEXO collaboration, and present the results of characterization 
of SiPMs from several vendors. In particular, we find that the 
photon detection efficiency at the peak of xenon scintillation light 
emission (175-178 nm) approaches the nEXO requirements for 
tested FBK and Hamamatsu devices. Additionally, the nEXO 
collaboration performed radio-assay of several grams of bare 
FBK devices using neutron activation analysis, indicating levels of 
40 k, and of the order of <0.15, (6.9 10“4 .1.3 10“2), 

and <0.11 mBq/kg, respectively. 

Index Terms —silicon photomultipliers, xenon detectors, pho¬ 
todetectors 


1. Introduction 

S EARCH for neutrinoless double-beta decay (DBD) is 
an active field of research with important implications 
for nuclear and particle physics. nEXO is a planned 5-ton 
liquid xenon detector aimed to probe the effective Majorana 
neutrino mass in the inverted mass hierarchy region. nEXO 
benefits from experience gained with the successful EXO-200 
detector Q, which utilizes ^200 kg of liquid xenon kept at a 
nominal pressure of :^147 kPa and a temperature of :^167 
K. To achieve maximum sensitivity, the nEXO detector is 
required to have an energy resolution of ^ ^1% at 2.458 
MeV (Q value of DBD in ^^^Xe). The experiment further 
needs to achieve unprecedented levels of background in order 
to be able to detect a weak effect. A detector containing the 
largest physically or fiscally feasible amount of the decaying 
substance maximizes the decay rate. 

In the EXO-200 experiment, the resolution is limited by the 
noise in the photodetector channels. EXO-200 uses Large Area 
Avalanche Photo-diodes (LAAPDs) operated at a gain of ^200 
and an equivalent noise of ^10 photoelectrons per readout 
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channel. Alternative light sensors that are gaining popularity in 
recent years are Silicon Photomultipliers (SiPMs). SiPMs are 
photodetectors under investigation for use in nEXO because 
of their high gain and expected low radioactive content. Due 
to large gains, on the order of 10^, SiPMs have the potential 
to yield much higher signal-to-noise ratios, to the point that 
it is not unreasonable to expect noise below 1 photoelectron 
equivalent. SiPM devices that are currently available need to 
exhibit the following characteristics: 

1) High photon detection efficiency (PDE) at 175-178 nm 
(peak of Xe scintillation light emission): The aim for overall 
efficiency for detecting the scintillation photons produced in 
liquid Xenon is 10%. Simulations show that this efficiency 
is achievable if the sensors’ PDE is at least 15%. Note that 
this includes loss of light due to refiection off the SiPM’s 
surface, which is large (^50%) when SiPMs are submerged in 
liquid xenon due to the large mismatch in indices of refraction 
at 175 nm between silicon (n ^ 0.8 + 1.9i Q) and xenon 
(n = 1.7 0). As a first step of the photodetector R&D 
we perform measurements in vacuum in order to more easily 
test different photodetectors. Once devices passing our first 
level tests (described in this paper) have been identified, their 
functionality in liquid xenon will be verified and refiectivity 
of the photodetectors will be measured as a function of angle. 

2) High Radiopurity: We demand that the content of the 
photodetectors has to be a sub-dominant contributor to the 
background budget, compared to other sources (e.g., the 
detector vessel). As a guide, we currently require that the 
photodetectors contain no more than ^ppt levels of and 
^^^Th, which corresponds to specific activities on the order of 
^10 yuBq/kg. In practice, any solution that involves packaging 
material added to the bare device, especially ceramic, is 
disadvantageous, compared to using bare detectors. 

3) Small values of photodetector nuisance parameters: 
Dark noise pulses start contributing to the dead time if their 
rate exceeds ^200 MHz when summing over all SiPMs. 
Assuming 4 m^ photo-coverage, the dark noise rate limit is 
50 MHz/m^ or 50 Hz/mm^. This dark noise rate is to be 
achieved at liquid xenon temperatures. Correlated avalanches 
in SiPMs are caused by cross-talk (CT) and after-pulses (AP). 
They start contributing to the energy resolution if the fraction 
of additional (correlated) avalanches created by one avalanche 
exceeds 0.2. 

4) Good photodetector electrical properties: Scintillation 
light is emitted isotropically, hence the spatial distribution 
of the light does not need to be sampled finely. Decreasing 
the number of readout channels also allows us to minimize 
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TABLE I 

Summary of SiPM parameters required by nEXO. Parameters in 

ITALIC ARE PREEERABLE, BUT NOT MANDATORY. 


Parameter 

Value 

Photo-detection efficiency at I75-I78 nm 
(without anti-reflective coating in gas/vacuum) 

>15% 

Radio-purity: 

232Th and 238u 

<10 pBq/kg 

Dark noise rate at -I00°C 

<50 Hz/mm^ 

After-pulse and cross-talk probability 

<20% 

Single photodetector active area 

>1 cm^ 

Gain fluctuations and electronic noise 

<0.1 p.e. 

Single photon timing resolution 

<10 ns 


the number of required feedthroughs. The size of single 
photodetector elements should be at least 1x1 cm^ (the 
current baseline design for nEXO assumes a unit readout cell 
of 5 cm^). Larger areas, up to 20x20 cm^, are preferable 
but are limited by the electronics noise that scales with the 
capacitance per channel. Therefore, low capacitance per unit 
area is desirable (e.g. <50 pF/mm^) j^. 

5) Several ns timing resolution: Good timing resolution 
is not essential for nEXO. Several ns single photon timing 
resolution would be sufficient. Nevertheless, single photon 
timing resolution might prove useful for other applications, 
by allowing better pulse shape discrimination and position 
reconstruction with light only using time of flight. 

6) Compatibility with liquid xenon: Very long electron life 
time (^10 ms) must be achieved in nEXO, due to the long drift 
distances in a large detector, which requires very low levels of 
electron absorbing contaminants. Therefore, the photodetectors 
and associated ancillary mechanical and electrical interfaces 
will have to be tested in liquid xenon together with the passive 
optical elements to ensure acceptable levels of out-gassing. 
Furthermore, operation of the complete light detection system 
must be demonstrated in liquid xenon. In practice, any solution 
that involves wavelength shifters, which are commonly used 
to enhance PDE in the VUV region, poses additional risks, 
compared to the operation of bare devices. 

Table |l| summarizes the required SiPM parameters. 

The challenges listed above are addressed in close collab¬ 
oration with SiPM manufacturer Fondazione Bruno Kessler 
(FBK). In addition, SiPMs from other manufacturers, such as 
Hamamatsu, KETEK, Excelitas, Zecotek, and RMD, are under 
investigation too. We currently operate two setups for testing 
and characterization of SiPMs, allowing measurements of key 
parameters at cryogenic temperatures. These two setups are 
described in the following sections along with current results 
obtained with the FBK, Hamamatsu, and KETEK samples. 
Details on the devices used in this study are listed in Table |II 
The pictures of representative SiPM samples are shown Fig [1 

H. Test setup for PDE measurement 

The test setup is based on the earlier version employed for 
testing of the APDs for the EXO-200 experiment The 


TABLE II 

SiPMs used in this study. Devices eor which PDE measurements 

HAVE BEEN PEREORMED ARE GIVEN IN ITALIC. 


Producer 

UV sensitive 

Bare device 

Comments 

FBK 

Yes 

Yes 

Developed in 2010 based 
on the original n+/p FBK 
SiPM technology j^. 
This type will be referred 
as FBK-2010 in the text 

Hamamatsu 

Yes 

No 

Developed in 2013 for 
the MEG experiment 0. 
This type will be referred 
as Hamamatsu-VUV in 
the text 

KETEK 

No 

No 


FBK 

Yes 

Yes 

Developed in 2013 based 
on the ’’RGB” FBK tech¬ 
nology 0. This type will 
be referred as FBK-2013 
in the text 

Hamamatsu 

No 

No 

Similar technology 

to MEG, no UV 
enhancements. This 

device will be referred to 
as Hamamatsu-VIS in the 

text 







Fig. I. Examples of SiPMs used in this study. Five 4x4 mm^ FBK- 
2010 SiPM devices connected in parallel (left). The FBK-2010 devices were 
developed in 2010 based on their original n+/p technology and feature thin 
silicon oxide window to facilitate UV light transmittance. Four 6x6 mm^ 
Hamamatsu-VUV devices connected in parallel (right) developed for the MEG 
experiment, optimized for UV detection. 

setup was upgraded and optimized for the determination of 
PDE at 175-178 nm for devices of different types with surface 
areas up to 2x2 cm^. The schematic diagram and a picture 
of the setup are shown in Fig. The main components of 
the set-up are the detector cage and the source assembly, 
both enclosed in a vacuum chamber. The detector cage is 
a cylindrical copper box that houses the photodetector to be 
tested. Different photodetectors can be mounted on the orifice 
(2x2 cm^) in the bottom plate of the box using custom inserts 
that ensure equal distance from the source to the detector 
surface and good thermal contact between the detector and 
the box. The box is cooled by liquid nitrogen boil-off gas 
passing through a copper tube that is soldered along the 
circumference of the box’s bottom plate. The temperature is 
regulated by two resistive heaters controlled with Omega 
PID controllers. The source assembly consists of a custom 
built xenon scintillation source, an opaque optical cavity 
(4.5 inch diameter), and a bandpass optical filter (2 inch 
diameter). The xenon scintillation source utilizes 83 nCi (at 
the moment of fabrication; roughly 60 fissions per second 
during measurements) of ^^^Cf electroplated onto a platinum 
surface by Eckert & Ziegler pQ| , which is then enclosed 
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Fig. 2. A cross-sectional view of the CAD model (left) and a photograph 
(right) of the Stanford test setup. 


in a miniature chamber filled with xenon gas at ^1 bar 
pressure. A similar source based on Am was also fabricated, 
but not used routinely. The use of a ^^^Cf-fission source 
has the advantage of order of magnitude larger scintillation 
light output, compared to an alpha source (we estimate the 
scintillation flux at the bottom plate of the detector box to be 
photons/mm^ per fission) . A picture of one of the two 
xenon scintillation sources is shown in Fig. Additionally, a 
blue LED is installed inside the vacuum chamber and can be 
used as an alternative light source of varying intensity. The 



Fig. 3. The custom-built xenon scintillation source. 


Optical GD bandpass filter. The filter’s transmission curve 
is centered at 180 nm (40 nm FWHM) to avoid bias due to 
potential re-emission of Tefion into longer wavelengths and 
contamination by the sub-dominant infra-red component of 
xenon scintillation. 

A. PDE measurement 

Signals from the photodetectors are routed from the detector 
cage via a vacuum feed-through using low capacitance fiexible 
cables and amplified by a Cremat GD charge sensitive pream¬ 
plifier and a Gaussian shaper. The signals are then digitized 
by a National Instruments digitizer. A Lab View based data 
acquisition program collects waveforms in either continuous 
or triggered mode, which are then reconstructed and analyzed 
using ROOT GD based custom software. Fig. shows an 
example spectrum obtained with the FBK-2010 SiPM device 
and the setup described above. 

To provide an absolute PDE reference we use a Hamamatsu 
R9875P photomultiplier tube (PMT) with calibrated quantum 
efficiency at 175 nm. The reference PMT has a synthetic silica 
window, an outer diameter of 15.9 mm (Cs-Te photocathode 
effective diameter 8 mm) and is mounted in the detector cage 
at the same position as other photodetectors. The PMT is 
operated at room temperature, biased to the recommended 
plateau voltage (1190 V), and its gain is calibrated with low 
intensity LED light pulses. The PDE of the SiPM devices is 
then extracted by comparing their responses (specifically, the 
^^^Cf spontaneous fission peak position, expressed in number 
of photoelectrons per unit area, which is typically in the range 
of 0.2 - 0.4 photoelectrons per mm^) to that of the reference 
PMT. 


10 ^ 

10 " 

10 

1 
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Fig. 4. Charge spectrum recorded with the FBK-2010 SiPM device exposed 
to xenon scintillation light excited by the (black line) and 

(red hatched line) sources. The dark noise distribution is shown as blue filled 
region. Peaks from alpha and fission decays are visible. undergoes both 

alpha and spontaneous fission decays, while ^^^Am only undergoes alpha 
decay. 



Peak integral (a.u.) 


optical cavity confines the xenon scintillation light emitted by 
the source, avoiding parasitic refiections off aluminum and 
Tefion surfaces inside the vacuum chamber. It is implemented 
as a Tefion (PTFE) sphere with two concentric openings for the 
source at the bottom and the detector at the top, approximately 
2 inch in diameter each. Fitted at the top sphere’s opening, 2-3 
mm in front of the detector’s surface, is the Pelham Research 


An important part of the measurement is the parasitic charge 
correction applied to the SiPM devices. The total amount of 
observed charge is biased towards larger values, predominantly 
due to the contribution of correlated avalanches (typically 
dominated by cross-talk, with addition of after-pulses). The 
correction is determined from separate estimates of the rates 
of cross-talk and after-pulsing. Note that it is also possible 
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for dark hits to contribute to the bias, but in practice this 
contribution is typically negligible. 

For the PDF measurement, several sources of errors are 
taken into account: 

1) Raw charge spectrum calibration: The gain calibration 
is applied for each value of the over-voltage (defined as the 
excess bias beyond the breakdown voltage of the device) to 
convert the raw charge into the number of photoelectrons per 
unit area. The over-voltage value itself depends on the value of 
the breakdown voltage, which is determined by extrapolating 
the gain-voltage dependence to zero gain using a linear fit. The 
relative SiPM gain at each over-voltage is extracted from the 
dark charge distribution. While SiPMs are known to typically 
have very good single photoelectron (s.p.e.) resolution, the 
devices measured here have large readout areas (^1 cm^, con¬ 
nected in parallel), which leads to increased peak widths due to 
a) an increased readout noise and b) an increased rate of pulse 
reconstruction errors related to pile-up. In practice, only single 
or double photoelectron peaks could be reliably resolved. The 
gain is then determined by fitting the dark pulse distribution 
with a function representing pedestal and each photoelectron 
peak by a Gaussian. Fig. shows the distribution of random 
pulses occurring in the absence of light (dark pulses) obtained 
with the FBK-2010 device. The fit error on the s.p.e. equivalent 
charge is typically on the order of 0.5-1.5%. This error also 



Fig. 5. Dark pulse charge distribution obtained with the FBK-2010 device. 
Data shown as dots. The red solid line is a two Gaussian fit to the pedestal 
and single photoelectron peaks. 

affects the determination of the breakdown voltage and is 
propagated by the linear fit of the gain-voltage dependence. 
Note that the breakdown voltage error (typically, ^0.1-0.2 
V) could be straightforwardly translated into an additional 
error on the gain at a given voltage, thus decoupling this 
uncertainty from the inherent dependence of the PDF on 
over-voltage. Given the linear dependence of the gain on 
the over-voltage, the ±0.2 V uncertainty on the breakdown 
voltage translates into ^4% gain error for the FBK-2010 
and ^0.6% for Hamamatsu-VUV devices, averaged over used 
over-voltages. Similarly, the imperfect temperature stability 
of the setup (typical RMS ~0.2-0.3°C) also contributes to 
the gain error via the over-voltage variation. To evaluate this 
contribution we determine the temperature dependence of the 
over-voltage. The FBK-2010 devices show approximately 8- 
10 V breakdown voltage variation between room temperature 
and -104°C. The rate of change of the breakdown voltage with 
temperature is itself a function of temperature, so we measure 


it close to the value of interest (Fig.[^ and find a slope of ^50 
mV/°C near -104°C. Conservatively assuming 1°C variation, 
the gain uncertainty contribution is ^1%. In total, a 5% error 



Fig. 6. SiPM gain as a function of bias voltage at -104°C (circles) and 
-114°C (squares), measured with the FBK-2010 device. The blue lines are 
linear fits (the best-fit parameters are shown on the plot). Statistical error on 
individual gain measurements is ~1%, which is comparable to the marker 
size. 

on the charge calibration is estimated. 

2) Solid angle: The setup is designed such that the de¬ 
tectors being tested and the reference PMT are located at 
the same distance from the source and their surface areas 
subtend a small fraction of one steradian (<0.03 sr). In one 
case (Hamamatsu-VUV) a non-standard detector holder had to 
be used, which resulted in an offset of approximately 3 mm 
for this detector relative to others. This is accounted for by an 
appropriate correction. In all cases, the position is known to 
1-2 mm that, given the distance between the source and the 
detectors, translates into ^3% uncertainty on the light fiux. 

3) Light source temperature dependence: Due to imperfect 
thermal isolation of the source assembly from the detector 
cage, the source temperature was observed to gradually de¬ 
crease over time after the detector cage is cooled down and 
maintained at the liquid xenon temperature (^104°C). The 
source temperature decreases ^4°C/hr from room temperature 
during the first two hours following the cool down. Coincident 
with the temperature decrease, the fission peak position was 
found to increase at the rate of ~4-5%/hr, while the SiPM 
gain remained constant, confirming that the SiPM remained in 
thermal equilibrium with the detector cage (Fig. |^. Conser¬ 
vatively assuming that measurements with different detectors 
are performed two hours apart after cool down, we assign a 
10% error to this effect. 

4) Fission peak equivalent charge: The fission peak equiva¬ 
lent charge is determined by a Gaussian fit (Fig. A typical 
fit error is <1%. While the resolution of a single device is too 
poor to observe any potential structures in the fission energy 
spectrum, the Gaussian provides a good fit, and its mean is 
much farther separated from the background than the alpha 
peak, which makes it a good variable to compare efficiency 
of different devices. 

5) PMT absolute reference: The total error budget is cur¬ 
rently dominated by the uncertainty on the absolute PDF 
reference, which in our case is provided by the reference 
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Fig. 7. Drift of the fission peak position over time (top) observed with 
a SiPM held at -104°C. The temperature of the scintillation source was 
found to decrease by ~4°C/hr from room temperature after the cool down 
of the detector cage. The SiPM gain (bottom) remained constant during the 
measurement, confirming that the SiPM was in thermal equilibrium with the 
detector cage. 


PMT described earlier. The PDE of a PMT can be expressed 
as a product of quantum and collection efficiencies. The 
quantum efficiency depends on the wavelength of the detected 
light and was measured by Hamamatsu to be 20.1% at 175 
nm for the PMT used (the calibration uncertainty is ^2%). 
The collection efficiency depends on the voltage difference 
between the photocathode and the first dynode. It typically 
varies from 80% to 100% Fig.4-12], but for the PMT 

type used in this study was reported by Hamamatsu to be 
70%. Hamamatsu did not provide a reference curve for the 
particular tube we used, so we approximate the collection 
efficiency by 70±10%. The central value is suggested by 
Hamamatsu, while the full variation from 03 Fig.4-12] is 
used as an uncertainty. The tube to tube variability is assumed 
negligible. An additional uncertainty arises from imperfectly 
known single photoelectron gain. The PMT gain is typically 
calibrated with a very low intensity FED source, such that 
the charge spectrum obtained by triggering on the FED’s 
driver pulse is dominated by pedestal and single photoelectrons 
events (the mean Poisson photoelectron rate is ^0.1). This 
way the contamination of multiple photoelectron events is 
negligible and the gain can be determined typically with 
sub-percent uncertainty. Due to the poor peak-to-valley ratio 
observed for the used PMT, we calibrated the gain using larger 
light intensity (with mean Poisson rates of ^2-3.5), such that 
the centroid of the resulting charge distribution is more clearly 
separated from the pedestal peak. A substantial remaining 
probability of the pedestal (i.e., zero photoelectron) events 
allows us to determine the mean Poisson rate, which provides 
the expected fractions of single-, double-, etc., photoelectron 
events. The fraction of the pedestal events is determined by a 
Gaussian fit, with the peak width and position fixed to the val¬ 
ues obtained in the dark conditions. No loss of pedestal event 


is expected, given the externally triggered threshold-less data 
acquisition. A 10% error on the Poisson rate is conservatively 
assumed. The single photoelectron gain is then determined 
from a multi-photoelectron fit (Fig. The stability of the fit 
was checked by repeating the calibration procedure at different 
PMT voltages, recovering an exponential dependence of the 
gain. The gain calibration was found to vary within 6% if the 
light intensity was changed such that the mean Poisson rate 
changed by a factor of ^2. We round the quadratic sum of 
the three contributions up to 19% overall uncertainty on the 
PMT reference. 



Fig. 8. Distribution of peak amplitudes (black dots) obtained with the 
reference PMT exposed to low intensity LED light. The distribution is fit by 
a function (blue solid line) representing the sum of the pedestal and several 
photoelectron peaks. The mean Poisson rate (~2) is determined by the fraction 
of pedestal events. 

6) Parasitic charge correction: Parasitic charge corrections 
for FBK-2010 and Hamamatsu-VUV devices are taken from 
published data ( for FBK-2010 and Q for Hamamatsu- 
VUV). In both cases, the correction takes into account the 
main sources of correlated noise: cross-talk and after-pulsing. 
The probability of either, as is pointed out in |[^, does not 
change with the incoming photon flux, or wavelength. FBK 
claims very small probability of correlated avalanches for the 
devices based on their original technology, even at high over¬ 
voltages (^5% at 5V over-voltage). For the Hamamatsu-VUV 
devices, on the other hand, the contribution is substantially 
larger (^30% at 2V over-voltage). Neither publication ex¬ 
plicitly states uncertainties on the measured correlated noise 
contributions to the total charge. In the case of Hamamatsu- 
VUV, the ±(2-4) abs.% spread observed in the correction 
values at each over-voltage for different devices is used as 
an estimate of the uncertainty. In the case of the FBK-2010 
devices, we performed a partial cross-check by evaluating 
the probability of cross-talk (which is the dominant source 
of parasitic charge) in situ from the ratio of single- to 
multi-photon signals in dark conditions. The dependence of 
the observed cross-talk probability on over-voltage (Fig. 
agrees with the parasitic charge dependence on over-voltage 
extracted from to within 1-2 abs.%. Given the values of the 
required corrections, these differences translate into ^6 and 
^2 rel.% error on the PDE for Hamamatsu-VUV and FBK- 
2010 devices, respectively. We conservatively use the larger 
value for all measurements. 
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Fig. 9. Probability of cross-talk, estimated from the ratio of single- to multi¬ 
photon events, as a function of over-voltage for a FBK-2010 device. 


TABLE III 

Systematic errors oe the PDE measurement. Items shown in 

ITALIC CORRESPONDS TO THE SOURCES OE ERROR CORRELATED AMONG 
DIEEERENT MEASURED DEVICES. 


Source of error 

Value, % 

PMT absolute reference 

19 

Light flux temperature dependence 

10 

Parasitic charge correction 

6 

SiPM charge calibration 

5 

Solid angle 

3 

Fission peak equivalent charge 

1 

Total 

-23 


Table [IIJ summarizes the systematic error budget. Items 
shown in italic correspond to errors correlated among different 
measured devices. The total error budget of the absolute 
efficiency measurement is estimated at ~20%. The relative 
performance of Hamamatsu-VUV and FBK-2010 SiPM de¬ 
vices is determined more accurately (^13%, dominated by 
the light flux temperature dependence and parasitic charge 
correction). 


III. Test setup for cross-talk and after-pulsing 

MEASUREMENT 

The dark noise, after-pulsing rates, and the cross-talk 
probability were measured in a test setup at TRIUMF. The 
setup was designed to rapidly measure the performances of 
photodetectors in conditions as close as possible to nEXO. The 
photodetectors were mounted on printed circuit boards that 
attach to a INSTEC HCP302 cooling chuck whose temperature 
is regulated by the combination of liquid nitrogen boil off 
gas running and a heater embedded in the chuck. The cooling 
chuck can easily reach -110°C even when surrounded by warm 
gas. The whole setup is enclosed in a light tight box filled 
with nitrogen gas whose purpose is to prevent condensation 
and to allow the propagation of 175 nm light. No light 
sources were used, however, for the data reported in this paper. 
All data reported in this paper were taken at -100°C cold 
chuck temperature. The photodetector temperature tracked the 
temperature of the chuck within a few degrees depending on 
the photodetector package. 

Due to several constraints, different devices were tested 
at TRIUMF and at Stanford (see Table |n]). The devices 
tested at TRIUMF were manufactured by KETEK, EBK and 


Hamamatsu. The KETEK devices are 3x3mm^, 50/im pitch 
SiPM. The EBK devices, EBK-2013, are Ixlmm^ second 
generation SiPMs. The Hamamatsu devices, Hamamatsu-VIS, 
are 3x3mm^, 50/im pitch non-VUV sensitive SiPMs (also 
called by Hamamatsu Multi-Pixel Photon Counters, MPPCs) 
whose nuisance parameters are expected to be very similar 
to the MPPCs developed for the MEG experiment. The per¬ 
formances of these devices were measured at about 5V over¬ 
voltage for the KETEK and EBK-2013 devices, and 2.5V for 
the Hamamatsu-VIS MPPC. 

Signals from the photodetectors were amplified by a fac¬ 
tor 100 using two Mini-Circuits MAR6-SM-I- amplifiers in 
series GD A 50 Ohm coaxial cable was used to connect the 
photodetector to the amplifier board. The MAR6-SM-F are very 
fast amplifiers that allow measuring the photodetector pulse 
shape with minimum distortions. Signals from the amplifier 
board were readout by a Lecroy Waverunner oscilloscope. 
Waveforms were recorded by triggering on signals above the 
electronics noise. The waveforms were 10 /is long with the 
sampling rate limited to 1 Giga-sample per second (GS/s) as 
a compromise between timing resolution and data processing 
speed. Examples of digitized waveforms are shown on Eig. [T^ 

Dark noise and after-pulse rates are measured by compiling 
the distribution of the time differences between the trigger 
pulse and the next pulse. At -100°C most of the waveforms 
have only one pulse and the timing difference is calculated 
using the time stamp of each waveform provided by the 
oscilloscope. Unfortunately, the oscilloscope has a dead time 
of 4 /is after each waveform that has to be accounted for when 
characterizing the timing distribution. After-pulses can occur 
within the same waveform as the trigger pulse. The waveforms 
are analyzed as follows: 1) a pulse finder identifies pulses 
whose shape and amplitude are inconsistent with electronics 
noise, 2) pulses are fitted using an analytical template repre¬ 
senting the pulse shape by minimizing an approximate chi- 
square that is calculated using the single bin noise (about 1 
mV), 3) if the chi-square is larger than a certain threshold the 
pulses are refitted by adding additional pulses until either the 
chi-square falls below the threshold or the chi-square does not 
improve significantly or the number of added pulses reaches 
a maximum of 5. The third procedure allows identifying 
partially overlapping pulses. The pulse template that fits all 
tested SiPMs is the convolution of a Gaussian distribution with 
one or two exponential functions.The pulse polarity is negative 
starting with a fast fall time and followed by a slower rise time. 
The fall time is calculated as the time it takes to for the pulse to 
go from 10% to 90%. The rising part of the pulse is described 
by exponential functions and parametrized by the exponential 
time constants. Two exponentials are typically needed when 
the parasitic capacitance across the quenching resistor is larger 
than a few femtofarad. 

The obtained distribution of the time differences between 
the trigger pulse and the next pulse is fit by a function 
to extract dark noise and after-pulse rates. The function is 
constructed assuming that both dark noise and after-pulsing 
affect the time difference distribution as described in (D- 
The fit function is expanded to account for several features 
that were either ignored or that are specific to the nEXO 
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Fig. 10. Examples of single avalanche pulse waveforms for a KETEK (left), FBK-2013 (middle), and Hamamatsu-VIS (right) SiPMs devices with the fit 
function superimposed. The function fits the prompt peak but it is truncated when displayed. 


setup. Recovery was ignored in |T^ , but it has a noticeable 
impact on the probability of generating and detecting after¬ 
pulsing. Indeed, the probability of triggering an avalanche of 
any kind depends on the over-voltage and our fit function 
assumes somewhat arbitrarily that it scales linearly with over¬ 
voltage, hence it follows the exponential recovery function. 
In the early part of the recovery, our pulse finder algorithm 
may miss pulses if their amplitudes are too small. We model 
this process by introducing a dead time that depends on the 
recovery time constant (150, 100, and 8 ns for the KETEK, 
EBK-2013 and Hamamatsu-VIS SiPMs, respectively). Our 
function ignores delayed cross-talk that do not suffer from 
recovery because their occurrence rate is relatively small. 
They could be included at the cost of yet another parameter 
accounting for delayed cross-talk probability. Our function 
also includes the oscilloscope deadtime between 8.5 and 12.5 
/iS. The complete function shown in ([T]) combines the functions 
shown in (|nl|, ([^, (|^, and The dark noise rate is handled 
by a single rate parameter, Rd. Each after-pulsing probability 
is parameterized by an exponential function with two free 
parameters - a time constant, and a probability, Pa^- 
The after-pulsing indexes are omitted in and © for 
simplicity. After-pulsing is approximated using a probability 
rather than using Poisson statistics for simplicity as well. 
This approximation is justified as long as all after-pulsing 
probabilities are low. is the recovery time constant, trd, 
tsds, and tsde are the recovery deadtime, oscilloscope start 
deadtime, and oscilloscope end deadtime respectively. 


i<nA 

Ptotal (t) = Pd (t) Pa^ (t) 
i =0 

i<nA 3<nA^3^i 

+PAt)Y.^PAA) n ^ 1 ) 

i=0 7=0 


PA{t) = 


[ 


0 , 

_ t 
1 — e '’'r 
PA 


e 




if t'^trd or tsds'^t'^tsde 

Otherwise 
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e ^A -j-e A —e pa— e ^ t>tsde 
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where A = (xr -h TAjIrr/TA 
Pnit) = 

Poit) = 



if t'^tsds 

0, 

if tsds<t<t. 

Pl^^-Rni't-'tsde+'tsds 

\ if t>tsde 


if t'^tsds 

g sds 

if tsds^t^tsde 

g (t tsde~^^sds) 

if t^tsde 


( 4 ) 

( 5 ) 


The cross-talk probability was inferred by measuring the 
fraction of single pixel avalanches in the trigger window. We 
define the cross-talk probability as the probability that a single 
pixel avalanche triggers one or more additional avalanches 
within 1 ns of the parent avalanche. The physical mechanism 
responsible for cross-talk is the absorption of photons pro¬ 
duced in the parent avalanche within the high field region of 
neighboring pixels. The same mechanism is known to lead 
to delayed avalanches, within a few ns up to several /as later 
due to the diffusion of the charge carriers produced by photon 
absorption in regions with zero or very small electric field. By 
construction, we only include prompt cross-talk and such a 
delayed cross-talk contributes to after-pulsing instead. Indeed, 
to measure cross-talk we calculate the probability pi that the 
trigger pulse is a single pixel avalanche and we define the 
cross-talk probability as 1 — pi. 


IV. Radio-purity measurements 

The nEXO group at the University of Alabama investigated 
the radioactivity content of EBK-2013 SiPMs by means of 
neutron activation analysis (NAA). Particular attention was 
paid to K, U and Th. The radio active decay series of 
^^^Th and are a background concern for nEXO. High 
energy gamma-rays, emitted during the decay of Th and U 
daughters, can penetrate the detectors interior and contribute 














TABLE IV 

Parameters oe the important 7 backgrounds eor nEXO fTT) . 


Parent Isotope 

half-life 

[ 1 - 107 ] 

Important decay 

Important 7 -ray 
[keV] 


1.248± 0.003 

4UfJ 

1460.822±0.006 

232Th 

14.0 ± 0.1 

20871 _>'^®®Pb 

2614.511 ± 0.10 

238u 

4.468 ± 0.003 

214Bi^214po 

2447.70 ± 0.3 


to the measured background. Potassium has a lower Q-value, 
so is not important for the neutrinoless DBD, but it can affect 
the accuracy of alternative physics searches. Table W shows 
the half-lives and 7 -ray energies for the important nEXO 
backgrounds. Direct gamma-ray counting of large samples 
allows us to measure Th, and U concentrations of a few 
hundred ppt. NAA allows for ppt and even sub-ppt sensitivities 
and requires smaller samples, than direct background counting, 
an important aspect for this study where kg-size samples are 
hard to obtain. However, the translation of NAA-determined 
Th and U concentrations into a background expectation value 
requires assumptions on the chain equilibrium, not required 
for direct background counting. Three samples obtained from 
FBK were activated at the MIT research reactor (MITR) in 
three separate beam times. In order to evaluate bulk and not 
surface properties of these devices, they were cleaned using 
the techniques described in fT^ . However, there were some 
variations from the technique described in and what was 
done here. The June 2013 activation sample was soaked for 
three days in 0.5M certified low Th/U nitric acid C3- During 
the preparation process the TiN backing to the FBK-2013 chips 
was lost. To counter the loss of the backing the acid soak was 
reduced to 12 h for the following activations. The integral 
thermal, epi-thermal, and fast fission neutron fluxes at MITR 
were determined by activating NIST certified fly ash l |20| (an 
element cocktail with known composition) and a sample of 
TiN. In earlier studies we established that sample port 2PH1 
(located close to the fuel element) has a high thermal neutron 
flux of about (4-5)-10^^ cm“^s“^ but also a relatively high 


^12 


-^s-i 


creation of significant Ti-related side activities that interfere 
with the Th and U analysis. These unwanted side activities are 
mainly due to the ^^Ti(n,p)'^^Sc and '^^Ti(n,p)^^Sc reactions. 
To suppress these backgrounds we activated the last two SiPM 
samples in MITRs sample port IPHl, showing a 230-fold 
reduced fast neutron fiux at the expense of a 6 -times lower 
thermal neutron fiux. We compensated for the lower thermal 
flux by means of longer exposure times, a viable solution 
given the relatively long live time of the Th and U activation 
products ^^^Pa and ^^^Np. For IPHl, the fast neutron fiux 
was found to be so low, that the fiy-ash could not provide 
a definitive measurement. To find the fast-neutron fiux a 
sample of TiN was also activated. After receiving the activated 
samples at the University of Alabama they were separated 
from the activation vials and then weighed to account for 
sample loss. Time and energy differential counting was then 
performed using two shielded Ge-detectors. The calibration of 
the gamma ray detection efficiency of the Ge detectors was 
performed using an Eckert & Ziegler |TQ| source solution. 


TABLE V 

Listing oe the samples analyzed in three activation campaigns 
AT MITR. We list the samples masses, MITR sample port, and 

COUNTING DELAY (ALTER END OE ACTIVATION). THE MAY 2014 
ACTIVATION WAS SPLIT INTO TWO TIME SEGMENTS AS THE SAMPLE WAS 
EJECTED EARLY. THIS IS TAKEN INTO ACCOUNT IN THE ANALYSIS. 


Activation 

Port 

Sample 
mass [g] 

Activation 
length [h] 

delay in 
counting [h] 

June 2013 

2PH1 

0.1705 

8.00 

103 

May 2014 

IPHl 

3.5519 

81.43 

127 

July 2014 

IPHl 

3.8906 

80.45 

33 


containing multiple gamma-ray emitters of known activity. 
The gamma-rays counted by the detector, in form of full 
absorptions peaks, were compared to the known activities to 
arrive at an energy dependent efficiency curve, described by 
a fifth degree polynomial fit. For the analysis of the samples, 
the code described in |T^ was used. Briefly, this allows the 
use of the full data set and produces a time dependent decay 
curve for each isotope. An example of the double differential 
energy and time analysis is shown in Fig. El The fitted 
sample activity at some reference time (No in Fig. [nj then 
summarizes all 7 -peak integrals evaluated at all times in form 
of a single value. For all sample fits the half lives are fixed to 
their tabulated values. The elemental concentrations were then 
inferred using properly averaged thermal/epi-thermal and fast 
neutron cross-sections HD- It was assumed that the relative 
isotopic abundances in the samples are given by the natural 
terrestrial abundance ratios. Table |V] shows the details of the 
three beam times performed with FBK-2013 chips. For the 
first two activation campaigns, sensitive analyses of K and U 
were hampered by the relatively long delay between end of 
irradiation and start of counting. 



Fig. 11. An example for the time dependence fit for ^^®Np, the blue dots are 
data and the red line is the best ht value. The half-life is fixed to its tabulated 
value (56.54 h), so only the activity at the end of the activation (67+/-1 Bq) 
is determined by the fit. 


V. Results 

A. Photon detection efficiency at 175-178 nm 

Fig. shows an example of the xenon scintillation spec¬ 
trum obtained with the Hamamatsu-VUV device in the Stan¬ 
ford setup. The fission fragment peak equivalent charge is 
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TABLE VI 

PDE RESULTS. Note an additional correction oe 15% has been 

ADDED IN THE CASE OE HAMAMATSU-VUV DEVICES, AS DESCRIBED IN 
THE TEXT. Relative uncertainty ~23%. 


Device 

Over-voltage, 

[V] 

Fission 

peak, 

[p.e./mm^] 

Parasitic 

charge 

correction, 

[%] 

PDE, 

[%] 

Hamamatsu-VUV 

1.5 

0.135 

14.4 

7.6 


1.7 

0.161 

19.2 

8.7 


1.9 

0.174 

24.0 

11.0 


2.2 

0.269 

31.3 

13.2 

FBK-2010 

3 

0.064 

2.2 

3.4 


4 

0.110 

4.4 

5.7 


5 

0.147 

6.7 

7.6 


6 

0.187 

8.9 

9.4 


determined with a Gaussian fit. The correction for parasitic 
charge is applied based on 0. 

Two additional corrections are needed in this particular case. 
First, the device’s packaging includes a 0.5 mm thick quartz 
window in front of the sensitive area. To compensate for 
the finite transmittance of the window, we apply a 10±1% 
correction Fig.5]. Second, due to specifics of the device’s 
packaging, the sensitive surface was positioned 3±1 mm 
farther away from the light source than other detectors and 
the reference PMT. This results in 5.Oil.5% correction for 
the smaller solid angle. 


Table [Vl| lists fission peak positions, correction values, and 
PDE numbers at different over-voltages for the measured 
devices. All measurements were performed at -104°C. 



Fig. 12. Charge spectrum from a Hamamatsu-VUV device obtained at 1.5 
V over-voltage. The horizontal axis represents the number of avalanches (true 
photoelectrons plus correlated cell discharges) per mm^. The data are shown 
as black dots. The blue solid line depicts a Gaussian fit to the fission fragment 
peak. 


Fig. shows the PDE as a function of over-voltage for 
the Hamamatsu-VUV and EBK-2010 devices. Note that while 
at similar over-voltages the EBK-2010 devices show several 
times lower PDE than Hamamatsu-VUV’s, they can tolerate 
much larger over-voltages due to substantially smaller rates of 
correlated avalanches, eventually reaching comparable PDE 
values. 



Fig. 13. PDE at 175-178 nm as a function of over-voltage for Hamamatsu- 
VUV (red squares) and FBK-2010 (blue dots) devices. Hatched areas represent 
the error bands. 


B. Photodetector nuisance parameters and timing analysis 


Eig. shows typical waveforms for the KETEK, EBK- 
2013, and Hamamatsu-VIS SiPMs along with the fitted func¬ 
tion. The parameters of the function that best reproduce each 
pulse shape are summarized in Table |VII| Hamamatsu-VIS 
MPPC exhibits a single rise time constant and a relatively 
slow fall time compare to the KETEK and EBK-2013 SiPMs. 
The KETEK and EBK-2013 SiPMs have similar pulse shapes 
with a very sharp peak first followed by a very slow decay 
time constant. Most of the charges produced in the avalanches 
is contained in the later part of the pulse shape. 

The cross-talk probabilities measured for all 3 devices are 
shown in Table |VII| The probabilities are remarkably similar 
for all 3 devices. They unfortunately all exceed 20%, which 
is our upper limit for the correlated avalanche probability that 
includes cross-talk and after-pulsing. The systematic errors on 
this measurement are small (on the order of a few percent) as 
it only requires proper identification of the single avalanche 
pulse which is very well defined for all devices 

As mentioned earlier, the distribution of the timing differ¬ 
ence between consecutive pulses is used to measure the dark 
noise and after-pulsing rates. The starting pulses are required 
to correspond to the oscilloscope trigger in order to properly 
account for the oscilloscope dead time. The starting pulses 
are also required to correspond to single pixel avalanches 
(i.e. triggers with cross-talk are excluded) in order to measure 
the after-pulsing rate generated by a single parent avalanche. 
On the other hand, the second pulse can have any amplitude 


(above the noise). Eig. shows the amplitude of the second 
pulse as a function of the time difference with the trigger 
pulse. The main band at ^40 mV corresponds to single 
pixel avalanches. Cross-talk yields pulses with amplitudes two, 
three, or more times larger. There is no correlation between 
time and amplitude except below 100 ns where a drop below 
the 40 mV value is clearly visible. After an avalanche the 
voltage across the diode indeed recovers with a time constant 
given by the product of the pixel capacitance and quenching 


















TABLE VII 

Operational, nuisance, and pulse shape parameters oe the SiPMs tested at -I04°C. 
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Ketek 

FBK-2013 

Hamamatsu-VIS 

Operating Voltage (V) 

29.50 

28.50 

59 

Breakdown Voltage (V) 

24.70 ± 0.05 

23.50 ± 0.05 

56.43 ± 0.05 

Pulse fall time (10 to 90%) (ns) 

0.9 

1.0 

3.5 

Pulse rise (exponential) time constant (ns) 

5.3 

2 

30.6 

Pulse 2nd fall time constant (ns) 

200 

300 

N/A 

Pulse 2nd fall time fraction 

0.14 

0.023 

N/A 

Crosstalk probability 

0.253 ±0.001 

0.274 ± 0.001 

0.249 ± 0.002 

Recovery time, Tr, (ns) 

750 ± 25 

360 ± 10 

21 ± 1 

Dark noise rate, Rd, (Hz/mm^) 

8.0 ± 0.2 

310 ± 1 

0.50 ± 0.01 

Total after-pulsing probability 

0.57 ± 0.02 

0.025 ± 0.001 

0.16 ± 0.01 

After-pulsing within 1 /xs 

0.18 ± 0.02 

0.016 ± 0.001 

0.033 ± 0.005 

After-pulsing time constant, (ns) / probability, Pai 

235 / 0.18 

179 / 0.01 

8.3 / 0.031 

2 nd After-pulsing time constant, ta 2 ^ (a^s) / probability, Pa 2 

2.5 / 0.04 

0.87 / 0.01 

8 / 0.01 

After-pulsing time constant, , (^s) / probability, PAg 

160 / 0.35 

8.2 / 0.006 

200 / 0.05 

4^^ After-pulsing time constant, ta^^, i/ds) / probability, Pa 4 

N/A/0 

N/A/0 

5000 / 0.07 



Time of next pulse (ns) 


Fig. 14. Statistical distribution of amplitude of the pulses following the 
trigger pulse as a function of the timing difference between this pulse and the 
trigger pulse for the Hamamatsu-VIS MPPC. The black line is the recovery 
function. 



Fig. 15. Probability that the next pulse occur between t and t+dt for 
the KETEK (red triangle), FBK-2013 (blue diamond) and Hamamatsu-VIS 
(black circle) SiPMs. The best ht functions are also shown. The gap in the 
distributions around 10^ is due to the oscilloscope deadtime. Data error bars 
omitted for clarity. 


resistance. The pulses that do not have lower amplitude even 
though they occurred within the recovery time scale must 
come from different pixels. It is likely that these pulses are 
delayed cross-talk, the seed charge carrier being created by 
photons from the parent avalanche subsequently diffusing to 
a neighboring pixel. This feature is also clearly visible for the 
FBK-2013 SiPM, while it is absent for KETEK. The recovery 
time constants are summarized in Table |VII| The recovery 


time is an order of magnitude faster for the Hamamatsu-VIS 
MPPC. 


The timing distributions are shown in Fig. \T5\ including the 
function that best fits the data. Each step on a log-log plot 
corresponds to a distinct exponential time dependence. The 
last step corresponds to distribution of dark pulses. It is clearly 
visible for all 3 devices, being largest for the FBK-2013 SiPMs 


and smallest for the Hamamamstu-VIS MPPC. The difference 
is striking because the FBK-2013 SiPMs are 9 times smaller in 
area than the KETEK and Hamamatsu-VIS devices. The dark 
noise rates per unit area are summarized in Table |VII| Both 
KETEK and Hamamatsu-VIS devices are below the required 
50Hz/mm^ specification. But the EBK-2013 SiPM dark noise 
rate is an order of magnitude too high (a similar value was 
found for the PBK-2010 device). After-pulsing rates were also 
extracted by fitting the distributions in Pig.[^ Very slow after¬ 
pulses are visible and to some extent they could be accounted 
for in the dark noise rate. In table |VIH we included the total 
after-pulsing rate and the after-pulsing rate within the first 10 
/iS that is relevant for our application. The after-pulsing rate 
is low for both PBK-2013 and Hamamatsu-VIS SiPMs but it 
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is high for KETEK devices, and exceeds the specification of 
20% for the combined rate of cross-talk and after-pulsing. 

C Radio-assay of FBK SiPMs 

Table IVIIII shows the concentrations of various elements in 
the EBK-2013 SiPMs as determined by NAA. 

It can be seen that the thorium concentration of all analyzed 
samples is no more than 0(10^ Pg/g)- However, there is 
variation of a factor of ten in the amount of thorium seen 
in the samples. This may be related to the amount of backing 
that was on the chip. Eurther analysis is being planned to 
understand chip to chip variations. Our sensitivity to uranium 
is currently limited by the amount of arsenic contained in the 
chips, which is used as a dopant. The interference is due to the 
large number of ^^As peaks and a similar half-life compared 
to 239Np. 

VI. Summary 

nEXO collaboration is investigating applicability of silicon 
photomultipliers for light detection in a large ultra-low back¬ 
ground liquid xenon experiment. Here we report results of 
characterization of the SiPM samples obtained from several 
vendors. In particular, for the first time, the PDE at 175 
nm is reported for EBK SiPMs. The device, based on a 
modified version of EBK’s original n+/p technology, shows an 
efficiency that approaches 10% at ^5V over-voltage, which is 
close to the value required by nEXO. We find that Hamamatsu 
devices optimized for light detection at 175 nm have PDE {in 
vacuo) are close to the required 15% at 2 V over-voltage, 
which is in a reasonable agreement with the measurements by 
the MEG collaboration Q. These results are encouraging for 
nEXO and other xenon experiments, which would benefit from 
using detectors capable of detecting scintillation light directly, 
without using wavelength shifters. 

Our measurements of the nuisance parameters show that 
cross-talk does not meet our specifications for the most recent 
available technologies of all three manufacturers. Hamamatsu 
MPPCs fulfill all other requirements. The tested EBK SiPMs 
show a too high dark noise rates, and KETEK SiPMs have 
also too high after-pulsing rates. Both Hamamatsu and EBK 
are working on improving the SiPM technology to address 
these issues. 

Einally, a sensitive radio-purity assessment was performed 
with a large sample of SiPMs, indicating levels of ^^K, 
^^^Th, and contamination at or below <0.15, (6.9-10“^ 
- 1.3-10“^), and <0.11 mBq/kg, respectively. In spite of the 
observed variability between different samples, these results 
are very promising. 


and advice. We thank Hamamatsu Photonics for providing 
samples of the MEG MPPCs and EBK for their continuous 
cooperation. 
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TABLE VIII 

Concentrations of the elemental constituents determined by NAA for three FBK-2013 SiPM samples: June 2013, May 2014, and 
July 2014. The value given is the best fit result with Gaussian statical errors and 10% systematic uncertainty associated with 

THE CALIBRATION OF THE COUNTING EFFICIENCY OF THE GERMANIUM DETECTORS. WHERE APPROPRIATE UPPER LIMITS DERIVED WITH THE METHOD 

OUTLINED IN ARE GIVEN. 



June 2013 irradiation 

May 2014 irradiation 

July 2014 irradiation 


Concentration 

Concentration 

Concentration 

Na[ng/g] 

69.2 ±3.9 ±7.0 

1.56 ±0.09 ±0.16 

7.11 ±0.39 ±0.72 

K[ng/g] 

-1.5±10.8±0.1 «17) 

-14.6±9.11±1.5 «4.8) 

3.81 ± 0.44 ±0.39 

Sc[pg/g] 

47.8 ±2.8 ±4.8 

25.7 ±1.4 ±2.6 

7.33 ±0.40 ±0.74 

Cr[ng/g] 

1.25 ±0.08 ±0.13 

3.75 ±0.21 ±0.38 

0.396 ±0.022 ±0.040 

Fe[ng/g] 

- 

34.4 ±1.9 ±3.5 

14.8 ±0.8 ±1.5 

Co[pg/g] 

686 ±37 ±69 

137 ±73 ±14 

43.4 ±2.4 ±44 

As[ng/g] 

288 ±14 ±29 

389 ±19 ±39 

590 ±29 ±60 

Br[pg/g] 

8.91 ±0.46 ±0.90 

96.6 ±4.8 ±9.7 

169±9±17 

Ti[Mg/g] 

(9.3±1.1±1.0)-10-3 

121 ±2 ±14 

38.3 ±1.6 ±3.9 

Sb[pg/g] 

146±10 ±15 

42.9±2.1 ±4.3 

16.1±0.8±1.7 

W[ng/g] 

-0.078 ± 0.013 ± 0.010 «0.02) 

1.31 ± 0.006 ± 0.14 

0.151 ±0.008 ±0.016 

Au[pg/g] 

27.7 ±1.3 ±2.8 

69.0 ±3.3 ±6.9 

19.0±1.0±1.9 

Ir[fg/g] 

10.4 ±33 ±1.1 «65) 

327 ±17 ±33 

16.1±1.8±1.7 

Th[pg/g] 

2.2 ±2.3 ±0.3 «6.0) 

3.15 ±0.23 ±0.32 

0.17±0.12±0.02 (o.nl^-JJ) 

U[pg/g] 

-8±12±1 «13) 

-7.18 ±4.21 ±0.72 «9.1) 

1.0±6.1±0.1 «10.9)‘ 








